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' A central quantity in mean-field magnetohydrodynamics is the mean electromotive force £, which in general depends on the mean 
\J^) ' magnetic field. It may however have a part independent of the mean magnetic field. Here we study an example of a rotating conducting 
body of turbulent fluid with non-zero cross-helicity, in which a contribution to £ proportional to the angular velocity occurs (Yoshizawa 
t | 1990). If the forcing is helical, it also leads to an a effect, and large-scale magnetic fields can be generated. For not too rapid rotation, 
i the field configuration is such that Yoshizawa's contribution to £ is considerably reduced compared to the case without a effect. In that 
case, large-scale flows are also found to be generated. 
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l— ~~ '■ 1 Introduction 
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■ Many studies of the large-scale magnetic fields in turbulent astrophysical bodies such as the Su n or the 



(N 



X 



i . uaiaxy are carriea out m tne irameworK or mean-neia electrody namics [see tne textDooKs by ivioriatt 
\ 19781 . Parker 19781 . Krause and Radler 198Cj . Zeldovich et al. 19831 ). It is based on the induction equation 
governing the magnetic field B, 



— = Vx(UxB- Wo J), (1) 



were U is the fluid velocity, J = V x B / fiQ the current density, r\ the magnetic diffusivity, and \xq the 
vacuum permeability. Both the magnetic field B and the velocity field U are considered sums of mean 
jjj ■ parts, B and U, defined as proper averages of the original fields, and fluctuations. The averages are 
assumed to satisfy the Reynolds averaging rules. The mean magnetic field B then obeys the mean-field 
induction equation 

ffD 

— = Vx(UxB + £- Wo J). (2) 



Here £ = u x b is the mean electromotive force resulting from the fluctuations of velocity and magnetic 
field, u = U — U and b = B — B. Generally, £ can be represented as a sum 

£ = £(0) + £(B) (3) 
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of a part £^ , which is independent of B, and a part £^ B ' vanishing with B. In many representations and 
applications of mean-field electrodynamics the part £^ of £ is ignored. Only the part £^ B \ which is of 
crucial importance for dynamo action, is taken into account. 

Here we focus our attention on the part £^ of £ . It may depend on non-magnetic quantities influencing 
the turbulence, in general also on U. If the magnitude of U is small, and if U varies only weakly in space 
and time, we may write 

F ( 0) = F (oo) + s . + (4) 

with <f| 00 ^ as well as Sy and T«j. being independent of U. Of course, the contribution to £^ 

can only be non-zero if the turbulence allows us to define a direction. For example, turbulence in a 
rotating body shows in general an anisotropy determined by the angular velocity fi, and might 
then be proportional to fi, say equal to c^fi. The S»j term in ([3]) can only be unequal to zero if the 
turbulence lacks Galilean invariance. In the case of isotropic turbulence it describes a contribution to £^> 
proportional to U, say equal to cjjU . Note that in forced turbulence Galilean invariance can be broken if, 
ind ependent of the flow, the forcin g is fixed in space and shows a finite correlation time (for an example 



sec 



Radler and Brandenburg] 2010l ). The term, if restricted to isotropic turbulence, corresponds to 



a contribution to £^ proportional to V x U, say equal to cw V x U. The coefficients cq and cw are, 
in contrast to cjj, pseud o scalars. The co ntributions cqCI and cwV x U to the mean electromotive force 
were first considered by Yoshizawal ( 1990l ). He found that both cq and cw are closely connected with the 



cross helicity u ■ b. In what follows the occurrence of the contributions c^fi and cvi/V x U to the mean 
electromotive force £ is called "Yoshizawa effe ct" . This effect has b een invoked to explain magnetic fields 



in accretion discs ( Yoshizawa and Yokoi 19931 ) and spiral galaxies ( Yokoi 



1996) . It has also been used t o 



explain the surprisingly high level of magnetic fields in young galaxies ([Brandenburg and Urpin 19981 ). 
because the amplification of the mean field by this effect is independent of any seed magnetic field owing 
to the field's linear as opposed to exponential growth with time. 
As for the part £^ B > of £, we recall here the traditional ansatz 

= aijBj + r] ijk B jtk . (5) 

It can be justified for cases in which B varies only weakly in space and time. In the simple case of isotropic 
turbulence it takes the form £( g ) = gB — r?tV x B, wh ich describes the a effect and the occurrence of a 
turbulent magnetic diffusivity (jKrause and E,adlerlll98nl ). 



In this paper we report on numerical simulations of magnetohydro dynamic turbulence in a rotating 
body, that is, under the influence of the Coriolis force. We present results for the mean electromotive force 
and discuss them in the light of the above remarks, focussing particular attention on the Yoshizawa effect. 



2 Model 

We consider forced magnetohydro dynamic turbulence of an electrically conducting, compressible, rotating 
fluid which is permeated by a magnetic field. An isothermal equation of state is used so that the pressure 
p and the mass density p are proportional to each other, p = pc^, with c s being a constant sound speed. 
The magnetic field B, the fluid velocity U and the mass density p are assumed to obey 

dA 

— = U x B - r]p J + f M , (6) 



Dl> -c*Vlnp-2nxU+~JxB + -V-2pvS + f K , (7) 



Dt p p 
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Dlnp 



Dt 



V • U . (8) 



Unless indicated otherwise, we exclude a homogeneous part of the magnetic field. A is the magnetic vector 
potential, V x A = B, and r/ again the magnetic diffusivity, D/Dt = d/dt + U ■ V is the advective 
time derivative, fi the angular velocity which defines the Coriolis force, Sy = ^(Uij + Ujj) — ^5ij V • U 
the trace-less rate of strain tensor, v the kinematic viscosity, while /m and /k define the magnetic and 
kinetic forcing specified below. The simultaneous magnetic and kinetic forcing is a simple way to generate 
non-zero cross helicity. We admit only small Mach numbers, that is, only weak compressibility effects. 

The equations ©-([8]) are numerically solved in a cubic domain with the edge length L assuming periodic 
boundary conditions. Then k\ = 2ir/L is the smallest possible wavenumber. We assume that fl is parallel 
to the positive z direction, that is, fi = (0, 0, ft) with Q > 0. 

The forcing is as sumed to be constant during each time step but changes from one step to the next 
(|Brandenburgll200lh . We define it until further notice by 



M 



N M Re{~f k (t) exp[ifc(t) • x + ty(t)]} , f K = N K Re{ik(t) x ~f k{t) exp[ifc(t) • x + , (9) 



with 



-^m = NMC s \/noPoc s /k{5t , N K = M K c s \/c s /ki5t , (10) 



where Nm and A/k are dimensionless amplitudes, po is the initial mass density, considered as uniform, kf 
the average forcing wavenumber and St the duration of the time step. Furthermore, 

~. _ fk(t) -iek(t) x f k(t) 

f m - j=p , (11) 

where k(t) = k(t)/\k(t)\ is the unit vector of k(t) and f k m is a statistically homogeneous isotropic non- 
helical random vector field (Haue Both forcing functions /m and /k are non-helical if s = 0, 



and maximally helical if |e| = 1. They are 5-correlated in time, so for a discrete time step 5t they must be 
proportional to St^ 1 ^ 2 . Owing to <5-cor relation, the fluid flow in our model is Galilean invariant. We put 

= k{t) x e(t) 

(t) ViHt) 2 - W) ■ e(t)) 2 )St ' 

where e(t) is a random unit vector not parallel to k(t). In this way we have V • /m = V • /k = 0. The 
wavevectors k are chosen such that their moduli k = \k\ lie in a band of width 5k around a mean forcing 
wavenumber kf, that is, kf — 5k < k < kf — 5k, and we choose 5k = k\. 

We describe our simulations using the magnetic Prandtl number Pr^-, the Coriolis number Co, the 
magnetic Reynolds number Rejvf, and the Lundquist number Lu, 

Pr M = zVr?, Co = 2n/u rms kf , Re M = u xms /r]kf , Lu = b rms /^fj, p r]kf , (13) 

with u rms and 6 rms being defined using averages over the full computational volume. While Pr^f and Co 
are input parameters, ReM and Lu are used for giving results. For our numerical simulations we use the 
Pencil CODE0, which is a high-order public domain code (sixth order in space and third order in time) 
for solving partial differential equations, including the hydromagnetic equations given above. 



: //pencil-code . googlecode . com/ 
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Figure 1. Non-helical case. Dependence of Rejif and Lu on Co. 

3 Results and Interpretation 

We have performed a series of simulations with Ptm = 1, JVk = 0.01, Mm = 0.005, kf = 5fci and varying 
Co. As initial conditions we used U = A = and p = po- 

We discuss the results here in terms of space averages taken over the full computational volume defined 
above. These averages, which we denote by angle brackets, are of course independent of space coordinates. 
We put first B = (B) + b and U = (U) +u. Using B = V x A and the periodicity of A, we have (B) = 0, 
that is, B = b. By contrast, (U) is not necessarily equal to zero. (B) = is however enough to justify 
(U ■ B) = (u ■ b) and (U x B) = (u x b). 

Within this framework the mean electromotive force discussed above and denoted there by £ is equal 
to (u x b). According to the ideas expressed in the Introduction, and recalling that volume averages of 
spatial derivatives of our periodic variables A or U vanish, we expect 



(uxb) = cn^t + cuiU) 



(14) 



with cq determined by the cross-helicity (u ■ b). Owing to Galilean invariance of the flow in our model cjj 
should vanish. In all simulations under the mentioned conditions (U) turned out very small. Even if the 
initial condition for U was changed and so larger \(U)\ were generated, no influence of (U) on (u x b) was 
observed. We conclude from this that indeed cjj = 0. 

Let us give further results first for non-helical forcing, e = 0. In this case we expect no a effect and 
see no reason for the generation of large-scale magnetic fields. Figure [T] gives ReM and Lu, considered as 
measures for i£rms and 6 rms , as functions of Co. Figure [2] shows that the cross helicity (u • b) and, if Co / 0, 
also the z component of the mean electromotive force (u x b) are non-zero. The moduli of the x and y 
components of (u x b) are negligible. According to Yoshizawa's result we expect (u x b) z = ^C( u ' b) Co 
with £ being a number of the order of unity. Figure [3] shows that (u x b) z /(u ■ b)Co is indeed around 0.5. 

Consider next the case of maximally helical forcing, e = 1. The simulations for this case have been 
carried out with a modified definition of /k- In Q and (fT0|) . ik(t) x f k M has been replaced by ffg(t)> 

and sj Cs/k[5t by a/ c^kf/di. Now an a effect is to be expected and, as a co nsequence, the gene ration of 
magnetic fields with scales comparable to that of the computational domain (|Brandenburd 200ll ). Indeed, 
as illustrated by Figure [H different typ es of large-scale magnetic fields with a dominant wavenumber k = k\ 
occur. Following Hubbard et all (120091 1. we call them "meso-scale fields" . As can be seen in the example of 
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Figure 2. Non-helical case. Normalized cross helicity (it • b) / u TrnB b TulB and z component of normalized mean electromotive force (u X 
b) /«rmsf>rms as functions of Co; the moduli of the x and y components of (it X b) /u rms b IIns are below 10 — 3 . 
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Figure 3. Non-helical case. Dependence of (it X 6} z /(it ■ b)Co on Co. 

Figure El these fields are to a good approximation of Beltrami shape. Three different types of such fields 
have been observed, 

B x = -Bo(0, sin k±x, cos k\x) , B Y = Bq(cos kiy, 0, sin k\y) , B z = f?o(sin k\z, cos k±z, 0) , (15) 

in general with common phase shifts of the components in the x, y and z directions. Bq was always of 
the order of several equipartition values B eq , defined by B eq = ^/J-opo n rms . For not too large Co all three 
types, B x , B Y and B z , turned out to be possible, but for Co exceeding a value of about unity only that 
of type B z occurs. This becomes understandable when considering that for the amplification of meso-scale 
fields of type B x and B , the products a vv a 77 , and a^ x a zz are important, while for B z it is a xx a yy , but 
\a zz \ is reduced by rotational quenching (Rudiger 1979 ) for large values of Co. 



7, 2011 1:33 Geophysical and Astrophysical Fluid Dynamics paper 

6 A. BRANDENBURG and K.-H. RADLER 




Figure 4. Helical case. Upper row: B y /B e<1 on the periphery of the computational domain, state with B x type field (left) and B z 
type field (right), Co = 0.37. Lower row: same as above, but U y /u Ims . 

Furthermore, meso-scale flows of type U x and U Y , defined analogously to (fT5l) . are also possible; see 
the lower panels of Figure [U Such flows have never been seen in the absence of cross helicity. How- 
ever, already for small Co it seems impossible to tolerate U z flows. This might be connected with the 
fact that the Coriolis force acting on a U z flow would produce a 90° phase-shifted flow proportional to 
(cos kiz, — sin kiz, 0). By comparison, the Coriolis force acting on a U or a U Y flow gives another one 
proportional to (sin k%x, 0, 0) or (0, cos k\y, 0), respectively, which does not directly interfere with U x or 
U Y . 

Both the cross helicity (u ■ b) and the mean electromotive force (u x b) are influenced by the presence of 
the meso-scale magnetic fields and meso-scale flows. Figure [6] shows the dependence of (u ■ b) and (u x b) z 
on the types of the meso-scale magnetic fields and on Co. Meso-scale magnetic fields of B x or B Y type 
together with meso-scale flows enhance the level of {u ■ b) /u Tms b rras , especially for small values of Co. With 
meso-scale magnetic fields of B z type (u ■ b)/u rms b rms is reduced relative to that in the non- helical case 
(Figure [2]), because 6 rms is enhanced by a factor of about 2. As Figure [3 demonstrates, (u x b) z /(u ■ b)Co 
depends now crucially on whether meso-scale fields of B x or B Y type or of B z type are present. In the 
first case the Yoshizawa effect is clearly reduced by the meso-scale fields; in the second case it is enhanced 
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Figure 5. Helical case. Profiles of B y (x)/ B C q and B z (x)/ B C q as well as their product in a state with B z field, Co = 0.2. Overbars 
denote yz averages. The dashed line gives the level of the x average of B y B z /B^ q , which is close to zero (here, fa -10" 3 ). 
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Figure 6. Helical case. Normalized cross helicity (it • b) /ti rms i>rms (upper lines) and z component of the normalized mean electromotive 
force (u X b) /u Tms 6 rms (lower lines) as functions of Co; the moduli of the x and y components are below 10 . Solid lines correspond to 
states with B x or B Y type fields, dashed lines to states with B z type fields. 



for small Co, but reduced for larger Co. 

The remarkable strength of the meso-scale fields can lead to strong magnetic quenching effects. As a 
first approach to the understanding of such effects the non-helical case has been studied with an imposed 
homogeneous magnetic field in the y or z directions, (0, Bq, 0) or (0, 0, Bq), respectively. Figure [8] shows as 
an example the dependence of (u x b) z /{u ■ b)Co at Co ~ 0.25 on Bo/B eq . It suggests that in the helical 
case the reduction of (u x b) z /(u ■ b)Co by B x or B Y fields, which possess a non-zero z component, is 
stronger than that by B z fields, which have no z components. 
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Figure 7. Helical case. Dependence of (u X b) z /{u • b)Co on Co. Solid lines correspond to states with B x or B Y type fields, dashed 
lines to states with B z type fields. 




Figure 8. Non-helical case with an imposed homogeneous magnetic field in y or in z direction, (0, Boi 0) (dashed line) or (0, 0, So) (solid 
line). Dependence of (u X b) z /{u ■ b)Co on Bo/B eq at Co ss 0.25, Rej^f ~ 10. 



4 Discussion 



The mean electromotive force in a turbulent fluid may have a part which is independent of the mean 
magnetic field and also independent of the mean flow. As an example for that we have studied forced hy- 
dromagnetic turbulence in a rotating body. In this case the Yoshizawa effect occurs, that is, a contribution 
cqQ to (tt x b). We have confirmed that cq is determined by the mean cross-helicity (u ■ b). We have 
also seen that, if an a effect is present, the Yoshizawa effect can to a large extend be compensated by the 
action of magnetic fields maintained by this a effect. 
In astrophysics, the occurrence of non-zero cross-helicity is not a very common phenomenon. We give 
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here a few examples in which the findings of this paper could be of interest. In the solar wind the systematic 
radial flow together with the Sun's large-scale magnetic field give rise to cross helicity of opposite sign in 
the two hemispheres. Although this primarily implies cross helicity associated with mean flow and mean 
magnetic field, it also results in cross helicity associated with the fluctuations. Together with the Sun's 
rotation, the latter should then produce a component of the mean electromotive force that is distinct 
from that related to the a effect. Note, however, that the cross-helicity associated with the fluctuations is 
directly a consequence of the cross helicity from the large-scale field. 

Another exa mple where small-sca le cross helicity can be generated is in a stratified layer with a vertical 
magnetic field (jRiidiger et q/.ll201ll ). Again, the sign of (u ■ b) is linked to the orientation of the large-scale 
field relative to the direction of gravity. 
Finally, cross helicity can be generated spontaneously and can then be of either sign, s uch a s 



in the Archontis dynamo ( Archontis! 2000T); for k i nemat ic simulations see lArchontis etall (|2003l ) 



Cameron and Gallowavl ( 2006 )7 ISur and Brandenburg! (2009) have analyzed this dynamo with respect to 



the Yoshizawa effect. In this example too, large-scale and small-scale fields are intimately related. This 
interrelation means that whenever we expect the £^ term to be present in an astrophysical system, there 
should also be a mean magnetic field. Such an effect that is odd in the mean magnetic field might therefore 
instead just as well be associated with an a effect. As it turns out, this is also the case in the present 
simulations, where a large-scale magnetic field has been produced. In the present case, we have gone a step 
further by including also kinetic helicity, in addition to just cross helicity. This produces an a effect and, 
as a consequence of this, a large-scale magnetic field. This field is particularly important when rotation is 
weak, because then the Yoshizawa effect is strongly quenched by this field. 
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